i=1
for(i in 1:length(dependentlist)){
dep=dependentlist[i]
print(dependentlist_eng[i])
load(paste("../data/modeldata/SVMorigmodeldatawithgeoandgeom_",dep,".RData",sep=""))
preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/SVM_geoandgeom/SVMwithgeoandgeom_fw_5fold_10p_",dep,"_allpreds",sep=""),kk=1:5,endround = 10,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
predict_radial_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,kappasum = F,tausum = F)
predictors = vector()
for (n in 1:ncol(preds)){
for (n2 in 1:nrow(preds)){
predictors <- c(predictors,as.character(preds[n2,n]))
}
}
print("###################### WITH PREDICTORS from the FW SELECTION ###################")
uniquepredictors <- unique(predictors)
predict_radial_full(modeldata=origmodeldata, dependent=dep, predictors=uniquepredictors)
print("######################WITH ALL PREDICTORS###################")
predict_radial_full(modeldata=origmodeldata, dependent=dep, predictors=unlist(paramsets[5]),printpreds = FALSE)
rm(origmodeldata,paramsets,paramsetnames,dependent)
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"
## [1] "Prediction error at end is: 0.646753246753247"
## [2] "Prediction error at end is: 0.637662337662338"
## [3] "Prediction error at end is: 0.675324675324675"
## [4] "Prediction error at end is: 0.657142857142857"
## [5] "Prediction error at end is: 0.648051948051948"
## [6] "Prediction error at end is: 0.657142857142857"
## [7] "Prediction error at end is: 0.648051948051948"
## [8] "Prediction error at end is: 0.620779220779221"
## [9] "Prediction error at end is: 0.620779220779221"
## [10] "Prediction error at end is: 0.610822510822511"
## [1] "10fold cv-error: 0.574074074074074"
## [1] "For predictors TRI_hr_ws25 AND geom_10m_fl4_L37"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 4 12 20 33 16
## 5 0 1 1 5 16
## [1] "Kappa overall = 0.172252533125487"
## [1] "Tau overall = 0.31712962962963"
## [1] "mean quality = 0.155384615384615"
## [1] "The quality is 0.155384615384615"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.518518518518519"
## [1] "For predictors TRI_hr_ws25 AND geom_10m_fl4_L37 AND TRI_hr_ws42 AND Slope_Height_hr AND profc_ws13_hr_hr AND longc_ws5_hr_hr AND Channel_Network_Base_Level_50m AND geom_hr_L50m_fl1_rplipshannon_UE_hr_10cells_hr AND Convexity_10m AND slope_DTM_50m_avg_ws9_50m AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND Standardized_Height_50m AND Normalized_Height_50m AND elev_ws15_10m AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_40cells_hr AND ProfileCurvature_10m AND minic_ws15_hr_hr AND fischerk_ws61_hr AND geom_10m_fl4_L9 AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_60cells_hr AND slope_DTM_50m_avg_ws7_50m AND Longitudinal_Curvature_hr AND dtm_hr_TEXTURE_r10 AND Texture_50m AND crosc_ws11_10m AND minic_ws5_hr_hr AND geom_dtm_10m_hyd_fl5_L70 AND geom_dtm_10m_hyd_fl5_L16 AND profc_ws5_hr_hr AND MRVBF_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr AND elev_ws11_10m AND crosc_ws3_10m AND crosc_ws3_hr_hr AND TRI_hr_ws18 AND fischerk_ws3_hr AND dtm_hr_TEXTURE_r20 AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_60cells_hr AND GeneralCurvature_10m AND Flow_Line_Curvature_hr AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_60cells_hr AND crosc_DTM_50m_avg_ws9_50m AND crosc_DTM_50m_avg_ws11_50m"
##
## preds 1 2 3 4 5
## 1 1 0 0 0 0
## 2 0 5 0 0 0
## 3 1 0 8 0 0
## 4 2 7 12 37 12
## 5 0 1 1 1 20
## [1] "Kappa overall = 0.501123595505618"
## [1] "Tau overall = 0.571759259259259"
## [1] "mean quality = 0.41816141604874"
## [1] "The quality is 0.41816141604874"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.569230450224967"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.564814814814815"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 0 1 0 0 0
## 3 0 0 0 0 0
## 4 4 11 20 36 17
## 5 0 1 1 2 15
## [1] "Kappa overall = 0.214443434212235"
## [1] "Tau overall = 0.351851851851852"
## [1] "mean quality = 0.178717948717949"
## [1] "The quality is 0.178717948717949"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"
## [1] "Prediction error at end is: 0.61991341991342"
## [2] "Prediction error at end is: 0.628571428571429"
## [3] "Prediction error at end is: 0.647619047619048"
## [4] "Prediction error at end is: 0.685714285714286"
## [5] "Prediction error at end is: 0.648484848484848"
## [6] "Prediction error at end is: 0.667099567099567"
## [7] "Prediction error at end is: 0.638095238095238"
## [8] "Prediction error at end is: 0.627705627705628"
## [9] "Prediction error at end is: 0.628138528138528"
## [10] "Prediction error at end is: 0.637662337662338"
## [1] "10fold cv-error: 0.481481481481481"
## [1] "For predictors geom_10m_fl8_L12 AND Longitudinal_Curvature_hr"
##
## preds 2 3 4 5
## 2 23 10 8 1
## 3 0 0 0 0
## 4 10 12 26 10
## 5 0 1 1 6
## [1] "Kappa overall = 0.294378698224852"
## [1] "Tau overall = 0.345679012345679"
## [1] "mean quality = 0.28653921687111"
## [1] "The quality is 0.28653921687111"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.481481481481481"
## [1] "For predictors geom_10m_fl8_L12 AND Longitudinal_Curvature_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr AND Catchment_Area_hr AND terraintexture_hr_ws33_tp25 AND terraintexture_hr_ws17_tp25 AND planc_DTM_50m_avg_ws3_50m AND geom_hr_L50m_fl1_rplipmps_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr AND slope_DTM_50m_avg_ws3_50m AND geom_10m_fl10_L27 AND vectorruggedness_hr_ws59 AND geom_hr_L3_fl1_rplipdominance_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_60cells_hr AND geom_10m_fl10_L80 AND StandardizedHeight_10m AND minic_ws5_10m AND Convergence_Index_hr AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L50m_fl1_rplipdominance_UE_hr_40cells_hr AND sagaTopographic_Wetness_Index_50m AND geom_10m_fl8_L9 AND Convexity_50m AND profc_ws11_hr_hr AND longc_ws5_hr_hr AND profc_ws3_10m AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_40cells_hr AND geom_10m_fl10_L70 AND Catchment_area_50m AND vectorruggedness_hr_ws35 AND VerticalDistancetoChannelNetwork_10m AND planc_DTM_50m_avg_ws11_50m AND MRRTF_hr AND geom_10m_fl3_L46 AND minic_ws3_hr_hr AND elev_ws11_10m AND minic_ws7_hr_hr AND Valley_Depth_hr AND geom_hr_L3_fl10_rplipmps_UE_hr_40cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND MassBalanceIndex_10m"
##
## preds 2 3 4 5
## 2 30 9 6 1
## 3 0 8 0 1
## 4 3 5 28 7
## 5 0 1 1 8
## [1] "Kappa overall = 0.55566311713456"
## [1] "Tau overall = 0.580246913580247"
## [1] "mean quality = 0.481657715717866"
## [1] "The quality is 0.481657715717866"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.58082023643375"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.592592592592593"
##
## preds 2 3 4 5
## 2 23 13 9 2
## 3 0 0 0 0
## 4 10 10 26 15
## 5 0 0 0 0
## [1] "Kappa overall = 0.201303584858361"
## [1] "Tau overall = 0.271604938271605"
## [1] "mean quality = 0.193734335839599"
## [1] "The quality is 0.193734335839599"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1a.3 - Habitat for soil organisms"
## [1] "Prediction error at end is: 0.511688311688312"
## [2] "Prediction error at end is: 0.482251082251082"
## [3] "Prediction error at end is: 0.453679653679654"
## [4] "Prediction error at end is: 0.425974025974026"
## [5] "Prediction error at end is: 0.389177489177489"
## [6] "Prediction error at end is: 0.398268398268398"
## [7] "Prediction error at end is: 0.380519480519481"
## [8] "Prediction error at end is: 0.416450216450216"
## [9] "Prediction error at end is: 0.361471861471861"
## [10] "Prediction error at end is: 0.389177489177489"
## [1] "10fold cv-error: 0.481481481481481"
## [1] "For predictors Channel_Network_Base_Level_hr AND Convexity_50m"
##
## preds 1 2 3 4
## 1 12 3 1 3
## 2 11 31 6 1
## 3 0 0 12 4
## 4 1 0 7 16
## [1] "Kappa overall = 0.532631578947368"
## [1] "Tau overall = 0.54320987654321"
## [1] "mean quality = 0.470812655086849"
## [1] "The quality is 0.470812655086849"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.569455556051712"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.296296296296296"
## [1] "For predictors Channel_Network_Base_Level_hr AND Convexity_50m AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND Slope_hr AND Texture_50m AND Maximal_Curvature_hr AND Valley_Depth_hr AND slope_DTM_50m_avg_ws5_50m AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND Normalized_Height_50m AND Slope_50m AND geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr AND terraintexture_hr_ws57_tp25 AND Normalized_Height_hr AND StandardizedHeight_10m AND fischerk_ws23_hr AND geom_dtm_10m_hyd_fl5_L6 AND elev_ws11_10m AND NormalizedHeight_10m AND TangentialCurvature_10m AND Catchment_slope_50m AND vectorruggedness_hr_ws7 AND Topographic_Wetness_Index_50m AND minic_ws7_hr_hr AND SGU AND Channel_Network_Base_Level_50m AND vectorruggedness_hr_ws3 AND minic_ws29_hr_hr AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_60cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr AND terraintexture_hr_ws45_tp25 AND maxic_ws7_10m AND geom_hr_L50m_fl1_rplipdominance_UE_hr_60cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr AND planc_DTM_50m_avg_ws9_50m AND slope_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws3_50m"
##
## preds 1 2 3 4
## 1 22 4 0 0
## 2 2 30 4 2
## 3 0 0 19 2
## 4 0 0 3 20
## [1] "Kappa overall = 0.787745664739884"
## [1] "Tau overall = 0.790123456790123"
## [1] "mean quality = 0.729828042328042"
## [1] "The quality is 0.729828042328042"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.800853249082389"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.564814814814815"
##
## preds 1 2 3 4
## 1 0 0 0 0
## 2 24 34 6 11
## 3 0 0 20 9
## 4 0 0 0 4
## [1] "Kappa overall = 0.346563407550823"
## [1] "Tau overall = 0.382716049382716"
## [1] "mean quality = 0.297857142857143"
## [1] "The quality is 0.297857142857143"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1a.4 - Habitat for crops"
## [1] "Prediction error at end is: 0.177056277056277"
## [2] "Prediction error at end is: 0.167099567099567"
## [3] "Prediction error at end is: 0.176190476190476"
## [4] "Prediction error at end is: 0.185281385281385"
## [5] "Prediction error at end is: 0.176190476190476"
## [6] "Prediction error at end is: 0.176190476190476"
## [7] "Prediction error at end is: 0.167099567099567"
## [8] "Prediction error at end is: 0.176190476190476"
## [9] "Prediction error at end is: 0.176190476190476"
## [10] "Prediction error at end is: 0.176190476190476"
## [1] "10fold cv-error: 0.166666666666667"
## [1] "For predictors slope_ws3_hr_hr AND slope_ws15_10m"
##
## preds 3 4 5
## 3 0 0 0
## 4 7 61 6
## 5 1 1 32
## [1] "Kappa overall = 0.719917012448133"
## [1] "Tau overall = 0.791666666666667"
## [1] "mean quality = 0.537777777777778"
## [1] "The quality is 0.537777777777778"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.138888888888889"
## [1] "For predictors slope_ws3_hr_hr AND slope_ws15_10m AND geom_hr_L3_fl10_rplipdominance_UE_hr_40cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_10cells_hr AND terraintexture_hr_ws53_tp5 AND terraintexture_hr_ws57_tp5 AND minic_ws11_10m AND Flow_Line_Curvature_hr AND SlopeHeight_10m AND Convexity_50m AND geom_10m_fl10_L6 AND TRI_hr_ws2 AND Slope_50m AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipedgedensity_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipshannon_UE_hr_60cells_hr AND geom_hr_L50m_fl10_rplipshannon_UE_hr_60cells_hr AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_60cells_hr AND slope_DTM_50m_avg_ws7_50m AND TRI_hr_ws57 AND Slope_hr AND slope_ws5_hr_hr AND geom_10m_fl1_L11 AND Convexity_10m AND geom_10m_fl10_L110 AND minic_DTM_50m_avg_ws3_50m AND profc_DTM_50m_avg_ws3_50m AND geom_10m_fl8_L3 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_40cells_hr AND TRI_hr_ws14 AND TRI_hr_ws9 AND LS_Factor_50m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_20cells_hr AND geom_10m_fl1_L4 AND TRI_hr_ws10 AND terraintexture_hr_ws37_tp25 AND terraintexture_hr_ws17_tp25 AND fischerk_ws7_hr AND TRI_hr_ws38 AND vectorstrength_hr_ws7_hr AND planc_DTM_50m_avg_ws3_50m AND geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_20cells_hr"
##
## preds 3 4 5
## 3 0 0 0
## 4 8 62 4
## 5 0 0 34
## [1] "Kappa overall = 0.775933609958506"
## [1] "Tau overall = 0.833333333333333"
## [1] "mean quality = 0.577524893314367"
## [1] "The quality is 0.577524893314367"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.185185185185185"
##
## preds 3 4 5
## 3 0 0 0
## 4 8 62 5
## 5 0 0 33
## [1] "Kappa overall = 0.75625"
## [1] "Tau overall = 0.819444444444444"
## [1] "mean quality = 0.565029239766082"
## [1] "The quality is 0.565029239766082"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.1 - Average precipitation retention capacity"
## [1] "Prediction error at end is: 0.565367965367965"
## [2] "Prediction error at end is: 0.593073593073593"
## [3] "Prediction error at end is: 0.593073593073593"
## [4] "Prediction error at end is: 0.62034632034632"
## [5] "Prediction error at end is: 0.620779220779221"
## [6] "Prediction error at end is: 0.648051948051948"
## [7] "Prediction error at end is: 0.648051948051948"
## [8] "Prediction error at end is: 0.648051948051948"
## [9] "Prediction error at end is: 0.657142857142857"
## [10] "Prediction error at end is: 0.638528138528139"
## [1] "10fold cv-error: 0.537037037037037"
## [1] "For predictors slope_ws15_10m AND GeneralCurvature_10m"
##
## preds 1 2 3 4 5
## 1 10 0 3 1 2
## 2 16 44 9 13 4
## 3 0 0 0 0 0
## 4 0 1 0 4 0
## 5 0 0 0 0 1
## [1] "Kappa overall = 0.273176761433869"
## [1] "Tau overall = 0.43287037037037"
## [1] "mean quality = 0.23432611701668"
## [1] "The quality is 0.23432611701668"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.518518518518519"
## [1] "For predictors slope_ws15_10m AND GeneralCurvature_10m AND TRI_hr_ws49 AND minic_ws15_10m AND ProfileCurvature_10m AND geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr AND terraintexture_hr_ws37_tp25 AND profc_DTM_50m_avg_ws5_50m AND fischerk_ws19_hr AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_20cells_hr AND crosc_ws5_10m AND Channel_Network_Base_Level_50m AND ValleyDepth_hr AND planc_ws13_hr_hr AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_10cells_hr AND crosc_ws11_10m AND elev_ws11_10m AND crosc_ws3_10m AND Plan_Curvature_hr AND CrossSectionalCurvature_hr AND Flow_Line_Curvature_50m AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr AND vectorruggedness_hr_ws57 AND geom_hr_L3_fl1_rpliprichness_UE_hr_60cells_hr AND profc_DTM_50m_avg_ws9_50m AND planc_ws7_hr_hr AND longc_ws3_10m AND Slope_50m AND slope_DTM_50m_avg_ws7_50m AND geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr AND LSFactor_10m AND RelativeSlopePosition_10m AND geom_10m_fl3_L50 AND LS_Factor_hr AND geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND maxic_ws15_hr_hr AND Slope_Height_50m AND geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND minic_ws3_hr_hr AND TRI_hr_ws42 AND TRI_hr_ws43"
##
## preds 1 2 3 4 5
## 1 18 0 3 0 3
## 2 8 45 9 12 3
## 3 0 0 0 0 0
## 4 0 0 0 6 0
## 5 0 0 0 0 1
## [1] "Kappa overall = 0.449865951742627"
## [1] "Tau overall = 0.560185185185185"
## [1] "mean quality = 0.324621212121212"
## [1] "The quality is 0.324621212121212"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.592592592592593"
##
## preds 1 2 3 4 5
## 1 3 0 0 1 1
## 2 23 45 12 17 6
## 3 0 0 0 0 0
## 4 0 0 0 0 0
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.0607334396289316"
## [1] "Tau overall = 0.305555555555555"
## [1] "mean quality = 0.10880721220527"
## [1] "The quality is 0.10880721220527"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.1 - Minimum precipitation retention capacity"
## [1] "Prediction error at end is: 0.805194805194805"
## [2] "Prediction error at end is: 0.769264069264069"
## [3] "Prediction error at end is: 0.776190476190476"
## [4] "Prediction error at end is: 0.767965367965368"
## [5] "Prediction error at end is: 0.777922077922078"
## [6] "Prediction error at end is: 0.758874458874459"
## [7] "Prediction error at end is: 0.758874458874459"
## [8] "Prediction error at end is: 0.768398268398268"
## [9] "Prediction error at end is: 0.741125541125541"
## [10] "Prediction error at end is: 0.693939393939394"
## [1] "10fold cv-error: 0.62037037037037"
## [1] "For predictors geom_hr_L50m_fl10_rplipdominance_UE_hr_40cells_hr AND vectorruggedness_hr_ws23"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 6 19 11 9 6
## 3 1 2 4 1 0
## 4 2 4 4 13 4
## 5 3 1 3 3 12
## [1] "Kappa overall = 0.27807486631016"
## [1] "Tau overall = 0.305555555555555"
## [1] "mean quality = 0.236286472148541"
## [1] "The quality is 0.236286472148541"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.648148148148148"
## [1] "For predictors geom_hr_L50m_fl10_rplipdominance_UE_hr_40cells_hr AND vectorruggedness_hr_ws23 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr AND CrossSectionalCurvature_hr AND geom_hr_L50m_fl10_rplipedgedensity_UE_hr_60cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND Total_Curvature_50m AND profc_ws11_10m AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND minic_ws11_10m AND planc_ws5_10m AND vectorruggedness_hr_ws33 AND minic_DTM_50m_avg_ws7_50m AND slope_ws5_hr_hr AND planc_ws29_hr_hr AND maxic_ws23_hr_hr AND planc_ws23_hr_hr AND crosc_ws11_10m AND geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr AND crosc_DTM_50m_avg_ws3_50m AND geom_hr_L50m_fl1_rplipmps_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_40cells_hr AND maxic_ws11_10m AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr AND Tangential_Curvature_50m AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND minic_ws15_10m AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_60cells_hr AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr AND Vertical_Distance_to_Channel_Network_50m AND profc_ws15_10m AND planc_ws11_10m AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_40cells_hr AND Relative_Slope_Position_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr AND Catchment_area_50m AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_40cells_hr AND vectorruggedness_hr_ws25 AND geom_hr_L50m_fl1_rplipdominance_UE_hr_40cells_hr AND minic_ws23_hr_hr AND longc_ws11_10m AND geom_hr_L3_fl10_rplipmps_UE_hr_10cells_hr"
##
## preds 1 2 3 4 5
## 1 5 0 0 0 0
## 2 1 21 2 3 2
## 3 1 1 16 1 0
## 4 4 4 2 21 3
## 5 1 0 2 1 17
## [1] "Kappa overall = 0.667180277349769"
## [1] "Tau overall = 0.675925925925926"
## [1] "mean quality = 0.573324283559578"
## [1] "The quality is 0.573324283559578"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.68251991023537"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.759259259259259"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 6 26 18 17 7
## 3 0 0 0 0 0
## 4 6 0 3 9 13
## 5 0 0 1 0 2
## [1] "Kappa overall = 0.13531799729364"
## [1] "Tau overall = 0.178240740740741"
## [1] "mean quality = 0.125161574618096"
## [1] "The quality is 0.125161574618096"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.2 - Retention capacity for heavy precipitation events"
## [1] "Prediction error at end is: 0.305627705627706"
## [2] "Prediction error at end is: 0.306493506493506"
## [3] "Prediction error at end is: 0.324242424242424"
## [4] "Prediction error at end is: 0.324242424242424"
## [5] "Prediction error at end is: 0.315151515151515"
## [6] "Prediction error at end is: 0.314718614718615"
## [7] "Prediction error at end is: 0.305627705627706"
## [8] "Prediction error at end is: 0.305627705627706"
## [9] "Prediction error at end is: 0.305627705627706"
## [10] "Prediction error at end is: 0.305627705627706"
## [1] "10fold cv-error: 0.268518518518518"
## [1] "For predictors crosc_ws7_10m AND maxic_ws11_10m"
##
## preds 1 2 3 4 5
## 1 75 10 6 3 8
## 2 0 1 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 5 0
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.264193792581378"
## [1] "Tau overall = 0.6875"
## [1] "mean quality = 0.29024064171123"
## [1] "The quality is 0.29024064171123"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.259259259259259"
## [1] "For predictors crosc_ws7_10m AND maxic_ws11_10m AND Tangential_Curvature_50m AND profc_ws19_hr_hr AND minic_DTM_50m_avg_ws7_50m AND Plan_Curvature_hr AND dtm_hr_CONVEX_r10 AND CatchmentArea_10m AND Catchment_Area2_50m AND Modified_Catchment_Area_50m AND maxic_ws15_10m AND Profile_Curvature_50m AND minic_ws5_10m AND Total_Curvature_50m AND longc_ws15_hr_hr AND geom_hr_L3_fl1_rplipedgedensity_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_40cells_hr AND Vertical_Distance_to_Channel_Network_50m AND terraintexture_hr_ws45_tp5 AND Minimal_Curvature_hr AND profc_ws11_10m AND Minimal_Curvature_50m AND longc_ws3_hr_hr AND profc_ws3_hr_hr AND crosc_ws11_10m AND Profile_Curvature_hr AND Convergence_Index_50m AND geom_10m_fl2_L8 AND Convexity_50m AND minic_ws11_hr_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_5cells_hr AND minic_ws13_hr_hr AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_10cells_hr AND planc_DTM_50m_avg_ws3_50m AND longc_ws11_10m AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND profc_ws5_10m AND geom_hr_L3_fl10_rplipshannon_UE_hr_20cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr"
##
## preds 1 2 3 4 5
## 1 75 11 6 1 6
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 7 0
## 5 0 0 0 0 2
## [1] "Kappa overall = 0.377969762419006"
## [1] "Tau overall = 0.722222222222222"
## [1] "mean quality = 0.376515151515152"
## [1] "The quality is 0.376515151515152"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.305555555555556"
##
## preds 1 2 3 4 5
## 1 75 11 6 8 8
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 0 0
## 5 0 0 0 0 0
## [1] "Kappa overall = 0"
## [1] "Tau overall = 0.618055555555555"
## [1] "mean quality = 0.138888888888889"
## [1] "The quality is 0.138888888888889"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.3 - groundwater reformation rate"
## [1] "Prediction error at end is: 0.637662337662338"
## [2] "Prediction error at end is: 0.647186147186147"
## [3] "Prediction error at end is: 0.665367965367965"
## [4] "Prediction error at end is: 0.646320346320346"
## [5] "Prediction error at end is: 0.60995670995671"
## [6] "Prediction error at end is: 0.609523809523809"
## [7] "Prediction error at end is: 0.637229437229437"
## [8] "Prediction error at end is: 0.638095238095238"
## [9] "Prediction error at end is: 0.647186147186147"
## [10] "Prediction error at end is: 0.656709956709957"
## [1] "10fold cv-error: 0.583333333333333"
## [1] "For predictors crosc_ws13_hr_hr AND crosc_ws3_10m"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 12 38 13 10 11
## 3 0 0 1 0 0
## 4 0 0 0 3 1
## 5 2 1 0 6 10
## [1] "Kappa overall = 0.23248730964467"
## [1] "Tau overall = 0.351851851851852"
## [1] "mean quality = 0.198213608023855"
## [1] "The quality is 0.198213608023855"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.574074074074074"
## [1] "For predictors crosc_ws13_hr_hr AND crosc_ws3_10m AND terraintexture_hr_ws25_tp25 AND dtm_hr_TEXTURE_r5 AND geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr AND geom_10m_fl8_L6 AND LongitudinalCurvature_10m AND geom_hr_L3_fl10_rplipmps_UE_hr_40cells_hr AND minic_ws3_hr_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr AND minic_ws19_hr_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipmps_UE_hr_10cells_hr AND Maximal_Curvature_50m AND crosc_ws5_10m AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_60cells_hr AND crosc_ws7_10m AND Channel_Network_Base_Level_hr AND planc_DTM_50m_avg_ws9_50m AND planc_ws3_hr_hr AND Flow_Line_Curvature_hr AND terraintexture_hr_ws41_tp5 AND CrossSectionalCurvature_50m AND minic_ws15_10m AND Longitudinal_Curvature_hr AND elev_ws15_10m AND terraintexture_hr_ws37_tp25 AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND Total_Curvature_50m AND crosc_DTM_50m_avg_ws5_50m AND geom_hr_L3_fl1_rplipmps_UE_hr_20cells_hr AND LS_Factor_hr AND vectorruggedness_hr_ws53 AND geom_hr_L3_fl1_rplipdominance_UE_hr_10cells_hr AND geom_hr_L3_fl10_rplipsimpson_UE_hr_40cells_hr AND ValleyDepth_10m AND Catchment_Area_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr AND profc_ws7_hr_hr"
##
## preds 1 2 3 4 5
## 1 6 0 0 0 0
## 2 8 39 7 4 7
## 3 0 0 4 0 0
## 4 0 0 2 15 0
## 5 0 0 1 0 15
## [1] "Kappa overall = 0.623286023574693"
## [1] "Tau overall = 0.664351851851852"
## [1] "mean quality = 0.536149068322981"
## [1] "The quality is 0.536149068322981"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.678740141625793"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.638888888888889"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 12 38 14 16 14
## 3 0 0 0 0 0
## 4 0 0 0 2 0
## 5 2 1 0 1 8
## [1] "Kappa overall = 0.158004158004158"
## [1] "Tau overall = 0.305555555555555"
## [1] "mean quality = 0.162591093117409"
## [1] "The quality is 0.162591093117409"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.4 - Potential for providing nutrients for plants"
## [1] "Prediction error at end is: 0.333766233766234"
## [2] "Prediction error at end is: 0.269264069264069"
## [3] "Prediction error at end is: 0.268398268398268"
## [4] "Prediction error at end is: 0.25021645021645"
## [5] "Prediction error at end is: 0.277922077922078"
## [6] "Prediction error at end is: 0.286147186147186"
## [7] "Prediction error at end is: 0.231168831168831"
## [8] "Prediction error at end is: 0.213419913419913"
## [9] "Prediction error at end is: 0.203896103896104"
## [10] "Prediction error at end is: 0.213419913419913"
## [1] "10fold cv-error: 0.287037037037037"
## [1] "For predictors minic_ws9_hr_hr AND Vertical_Distance_to_Channel_Network_50m"
##
## preds 1 3 5
## 1 58 5 8
## 3 0 0 0
## 5 7 6 24
## [1] "Kappa overall = 0.521227621483376"
## [1] "Tau overall = 0.638888888888889"
## [1] "mean quality = 0.425641025641026"
## [1] "The quality is 0.425641025641026"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.185185185185185"
## [1] "For predictors minic_ws9_hr_hr AND Vertical_Distance_to_Channel_Network_50m AND ValleyDepth_10m AND geom_10m_fl10_L13 AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl10_rplipsimpson_UE_hr_40cells_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_5cells_hr AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND planc_DTM_50m_avg_ws7_50m AND minic_ws5_hr_hr AND crosc_DTM_50m_avg_ws5_50m AND maxic_ws3_hr_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_40cells_hr AND planc_ws5_hr_hr AND geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_20cells_hr AND Normalized_Height_50m AND geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr AND SGU AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_40cells_hr AND ValleyDepth_hr AND crosc_ws29_hr_hr AND Total_Curvature_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_20cells_hr AND Mid_Slope_Positon_50m AND Topographic_Wetness_Index_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND terraintexture_hr_ws25_tp25 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_5cells_hr AND TRI_hr_ws1 AND Modified_Catchment_Area_50m AND geom_hr_L50m_fl1_rplipshannon_UE_hr_60cells_hr AND fischerk_ws23_hr AND dtm_hr_TEXTURE_r5 AND Profile_Curvature_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_5cells_hr AND crosc_ws3_hr_hr"
##
## preds 1 3 5
## 1 64 7 1
## 3 0 0 0
## 5 1 4 31
## [1] "Kappa overall = 0.759259259259259"
## [1] "Tau overall = 0.819444444444444"
## [1] "mean quality = 0.571516722201654"
## [1] "The quality is 0.571516722201654"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.416666666666667"
##
## preds 1 3 5
## 1 65 11 23
## 3 0 0 0
## 5 0 0 9
## [1] "Kappa overall = 0.256830601092896"
## [1] "Tau overall = 0.527777777777778"
## [1] "mean quality = 0.312605218855219"
## [1] "The quality is 0.312605218855219"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.5 - Potential as a CO2 sink"
## [1] "Prediction error at end is: 0.325108225108225"
## [2] "Prediction error at end is: 0.335064935064935"
## [3] "Prediction error at end is: 0.287878787878788"
## [4] "Prediction error at end is: 0.278787878787879"
## [5] "Prediction error at end is: 0.296969696969697"
## [6] "Prediction error at end is: 0.278354978354978"
## [7] "Prediction error at end is: 0.269264069264069"
## [8] "Prediction error at end is: 0.306060606060606"
## [9] "Prediction error at end is: 0.325108225108225"
## [10] "Prediction error at end is: 0.297402597402597"
## [1] "10fold cv-error: 0.305555555555556"
## [1] "For predictors Channel_Network_Base_Level_hr AND crosc_ws11_hr_hr"
##
## preds 1 2 3 4 5
## 1 44 0 6 1 0
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 7 1 12 32 5
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.518796992481203"
## [1] "Tau overall = 0.62962962962963"
## [1] "mean quality = 0.262068965517241"
## [1] "The quality is 0.262068965517241"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.268518518518518"
## [1] "For predictors Channel_Network_Base_Level_hr AND crosc_ws11_hr_hr AND Topographic_Wetness_Index_50m AND Catchment_area_50m AND vectorruggedness_hr_ws11 AND terraintexture_hr_ws13_tp25 AND ValleyDepth_hr AND Catchment_Area_hr AND crosc_DTM_50m_avg_ws9_50m AND Catchment_Area2_50m AND elev_ws11_10m AND maxic_ws3_hr_hr AND maxic_ws15_hr_hr AND MaximalCurvature_10m AND ValleyDepth_10m AND Mid_Slope_Positon_50m AND vectorruggedness_hr_ws5 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr AND slope_DTM_50m_avg_ws5_50m AND planc_DTM_50m_avg_ws9_50m AND terraintexture_hr_ws57_tp25 AND Vertical_Distance_to_Channel_Network_hr AND profc_DTM_50m_avg_ws7_50m AND planc_ws15_hr_hr AND planc_ws19_hr_hr AND vectorruggedness_hr_ws3 AND Catchment_slope_50m AND LSFactor_10m AND slope_DTM_50m_avg_ws7_50m AND ChannelNetworkBaseLevel_10m AND Valley_Depth_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND longc_DTM_50m_avg_ws5_50m AND geom_hr_L3_fl1_rpliprichness_UE_hr_60cells_hr AND minic_ws15_10m AND crosc_ws13_hr_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND maxic_DTM_50m_avg_ws7_50m AND longc_ws3_10m"
##
## preds 1 2 3 4 5
## 1 49 0 1 0 1
## 2 0 1 0 0 0
## 3 1 0 14 2 1
## 4 1 0 3 31 3
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.812399786210583"
## [1] "Tau overall = 0.849537037037037"
## [1] "mean quality = 0.667178387650086"
## [1] "The quality is 0.667178387650086"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.416666666666667"
##
## preds 1 2 3 4 5
## 1 47 0 6 7 2
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 4 1 12 26 3
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.458762886597938"
## [1] "Tau overall = 0.594907407407407"
## [1] "mean quality = 0.240537449971412"
## [1] "The quality is 0.240537449971412"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.1 - Potential for retention of heavy metals"
## [1] "Prediction error at end is: 0.417748917748918"
## [2] "Prediction error at end is: 0.380519480519481"
## [3] "Prediction error at end is: 0.38008658008658"
## [4] "Prediction error at end is: 0.38008658008658"
## [5] "Prediction error at end is: 0.388311688311688"
## [6] "Prediction error at end is: 0.351515151515152"
## [7] "Prediction error at end is: 0.37012987012987"
## [8] "Prediction error at end is: 0.351082251082251"
## [9] "Prediction error at end is: 0.36017316017316"
## [10] "Prediction error at end is: 0.36969696969697"
## [1] "10fold cv-error: 0.324074074074074"
## [1] "For predictors vectorruggedness_hr_ws47 AND dtm_hr_CONVEX_r30"
##
## preds 1 2 3 4 5
## 1 38 4 5 2 6
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 2 0
## 5 4 3 5 3 36
## [1] "Kappa overall = 0.519866629619339"
## [1] "Tau overall = 0.62962962962963"
## [1] "mean quality = 0.312272205938575"
## [1] "The quality is 0.312272205938575"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.287037037037037"
## [1] "For predictors vectorruggedness_hr_ws47 AND dtm_hr_CONVEX_r30 AND TWI_10m AND Topographic_Wetness_Index_50m AND profc_ws3_hr_hr AND minic_DTM_50m_avg_ws9_50m AND dtm_hr_TEXTURE_r20 AND geom_hr_L50m_fl10_rplipmps_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_60cells_hr AND Catchment_area_50m AND vectorruggedness_hr_ws31 AND minic_ws5_hr_hr AND terraintexture_hr_ws45_tp25 AND Standardized_Height_50m AND planc_DTM_50m_avg_ws11_50m AND Tangential_Curvature_50m AND NormalizedHeight_10m AND vectorruggedness_hr_ws55 AND vectorruggedness_hr_ws45 AND minic_ws29_hr_hr AND TRI_hr_ws11 AND geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl1_rplipshannon_UE_hr_60cells_hr AND vectorruggedness_hr_ws49 AND geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr AND longc_ws19_hr_hr AND minic_ws15_10m AND maxic_ws5_10m AND TRI_hr_ws20 AND Texture_10m AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND dtm_hr_TEXTURE_r10 AND Modified_Catchment_Area_50m AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_40cells_hr AND dtm_hr_TEXTURE_r5 AND geom_hr_L3_fl1_rplipedgedensity_UE_hr_40cells_hr AND Channel_Network_Base_Level_hr AND slope_DTM_50m_avg_ws5_50m AND longc_ws11_hr_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr AND longc_ws3_10m AND vectorruggedness_hr_ws9 AND geom_hr_L3_fl10_rplipdominance_UE_hr_40cells_hr AND MaximalCurvature_10m AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND Longitudinal_Curvature_hr"
##
## preds 1 2 3 4 5
## 1 42 6 5 3 3
## 2 0 0 0 0 0
## 3 0 0 3 0 0
## 4 0 0 0 2 0
## 5 0 1 2 2 39
## [1] "Kappa overall = 0.674252810529202"
## [1] "Tau overall = 0.74537037037037"
## [1] "mean quality = 0.4254731853073"
## [1] "The quality is 0.4254731853073"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.37037037037037"
##
## preds 1 2 3 4 5
## 1 38 4 6 3 7
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 0 0
## 5 4 3 4 4 35
## [1] "Kappa overall = 0.46969696969697"
## [1] "Tau overall = 0.594907407407407"
## [1] "mean quality = 0.24538766270515"
## [1] "The quality is 0.24538766270515"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.2 - Potential for transforming organic contaminants"
## [1] "Prediction error at end is: 0.51991341991342"
## [2] "Prediction error at end is: 0.527705627705628"
## [3] "Prediction error at end is: 0.593506493506494"
## [4] "Prediction error at end is: 0.602164502164502"
## [5] "Prediction error at end is: 0.565367965367965"
## [6] "Prediction error at end is: 0.583549783549784"
## [7] "Prediction error at end is: 0.555844155844156"
## [8] "Prediction error at end is: 0.547186147186147"
## [9] "Prediction error at end is: 0.564935064935065"
## [10] "Prediction error at end is: 0.518614718614719"
## [1] "10fold cv-error: 0.5"
## [1] "For predictors Channel_Network_Base_Level_hr AND crosc_ws13_hr_hr"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 0 4 1 0 3
## 3 8 13 24 2 5
## 4 0 0 0 1 0
## 5 2 5 3 6 31
## [1] "Kappa overall = 0.367032967032967"
## [1] "Tau overall = 0.444444444444444"
## [1] "mean quality = 0.251433011433011"
## [1] "The quality is 0.251433011433011"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.453703703703704"
## [1] "For predictors Channel_Network_Base_Level_hr AND crosc_ws13_hr_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_40cells_hr AND Relative_Slope_Position_hr AND slope_DTM_50m_avg_ws3_50m AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L3_fl10_rplipdominance_UE_hr_10cells_hr AND Slope_50m AND Catchment_area_50m AND geom_hr_L50m_fl10_rplipmps_UE_hr_10cells_hr AND Topographic_Wetness_Index_50m AND maxic_DTM_50m_avg_ws9_50m AND dtm_hr_CONVEX_r20 AND longc_ws11_10m AND General_Curvature_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr AND geom_10m_fl3_L30 AND elev_ws11_10m AND TRI_hr_ws57 AND vectorruggedness_hr_ws55 AND Profile_Curvature_hr AND Standardized_Height_hr AND Texture_10m AND geom_10m_fl10_L23 AND Standardized_Height_50m AND ChannelNetworkBaseLevel_10m AND geom_hr_L50m_fl1_rplipdominance_UE_hr_40cells_hr AND minic_DTM_50m_avg_ws3_50m AND Mid_Slope_Positon_50m AND terraintexture_hr_ws45_tp25 AND LS_Factor_50m AND longc_ws3_hr_hr AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_5cells_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr AND Convexity_50m AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr AND TRI_hr_ws36 AND geom_10m_fl3_L34"
##
## preds 1 2 3 4 5
## 1 3 0 0 0 0
## 2 3 14 0 0 0
## 3 4 6 27 2 4
## 4 0 0 0 2 0
## 5 0 2 1 5 35
## [1] "Kappa overall = 0.651237890204521"
## [1] "Tau overall = 0.6875"
## [1] "mean quality = 0.488107887384483"
## [1] "The quality is 0.488107887384483"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.632451537505577"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.574074074074074"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 0 0 0 0 0
## 3 6 15 18 2 3
## 4 0 0 0 0 0
## 5 4 7 10 7 36
## [1] "Kappa overall = 0.265120967741935"
## [1] "Tau overall = 0.375"
## [1] "mean quality = 0.174129353233831"
## [1] "The quality is 0.174129353233831"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"
## [1] "Prediction error at end is: 0.0367965367965368"
## [2] "Prediction error at end is: 0.0367965367965368"
## [3] "Prediction error at end is: 0.0458874458874459"
## [4] "Prediction error at end is: 0.0367965367965368"
## [5] "Prediction error at end is: 0.0367965367965368"
## [6] "Prediction error at end is: 0.0367965367965368"
## [7] "Prediction error at end is: 0.0458874458874459"
## [8] "Prediction error at end is: 0.0458874458874459"
## [9] "Prediction error at end is: 0.0367965367965368"
## [10] "Prediction error at end is: 0.0367965367965368"
## [1] "10fold cv-error: 0.0185185185185185"
## [1] "For predictors PlanCurvature_10m AND planc_ws11_hr_hr"
##
## preds 4 5
## 4 2 0
## 5 2 104
## [1] "Kappa overall = 0.658227848101266"
## [1] "Tau overall = 0.962962962962963"
## [1] "mean quality = 0.740566037735849"
## [1] "The quality is 0.740566037735849"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.518568491248852"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.0370370370370371"
## [1] "For predictors PlanCurvature_10m AND planc_ws11_hr_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_60cells_hr AND fischerk_ws7_hr AND vectorstrength_hr_ws7_hr AND Convexity_10m AND Catchment_Area2_50m AND planc_ws9_hr_hr AND planc_ws13_hr_hr AND fischerk_ws3_hr AND crosc_ws11_10m AND crosc_ws15_10m AND crosc_ws3_10m AND crosc_ws5_10m AND crosc_ws7_10m AND longc_ws11_10m AND longc_ws15_10m AND longc_ws3_10m AND longc_ws5_10m AND longc_ws7_10m AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplippatchnum_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplippatchnum_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND planc_ws5_hr_hr"
##
## preds 4 5
## 4 3 0
## 5 1 104
## [1] "Kappa overall = 0.852459016393442"
## [1] "Tau overall = 0.981481481481481"
## [1] "mean quality = 0.870238095238095"
## [1] "The quality is 0.870238095238095"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.712718059941681"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.0370370370370371"
##
## preds 4 5
## 4 0 0
## 5 4 104
## [1] "Kappa overall = 0"
## [1] "Tau overall = 0.925925925925926"
## [1] "mean quality = 0.481481481481482"
## [1] "The quality is 0.481481481481482"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.4 - Potential for retention of water-soluble contaminants"
## [1] "Prediction error at end is: 0.592207792207792"
## [2] "Prediction error at end is: 0.555411255411255"
## [3] "Prediction error at end is: 0.564502164502165"
## [4] "Prediction error at end is: 0.601731601731602"
## [5] "Prediction error at end is: 0.601731601731602"
## [6] "Prediction error at end is: 0.592640692640693"
## [7] "Prediction error at end is: 0.573593073593074"
## [8] "Prediction error at end is: 0.564502164502165"
## [9] "Prediction error at end is: 0.564502164502165"
## [10] "Prediction error at end is: 0.573593073593074"
## [1] "10fold cv-error: 0.453703703703704"
## [1] "For predictors TRI_hr_ws29 AND fischerk_ws55_hr"
##
## preds 1 2 3 4 5
## 1 45 10 7 11 1
## 2 0 0 0 0 0
## 3 2 0 7 3 3
## 4 0 4 3 8 3
## 5 0 0 0 0 1
## [1] "Kappa overall = 0.325627740135512"
## [1] "Tau overall = 0.456018518518518"
## [1] "mean quality = 0.249421052631579"
## [1] "The quality is 0.249421052631579"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.416666666666667"
## [1] "For predictors TRI_hr_ws29 AND fischerk_ws55_hr AND geom_hr_L3_fl1_rplipshannon_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr AND slope_DTM_50m_avg_ws5_50m AND Normalized_Height_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND Catchment_area_50m AND Profile_Curvature_hr AND FlowLineCurvature_10m AND slope_ws11_hr_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_10cells_hr AND RelativeSlopePosition_10m AND crosc_ws7_hr_hr AND TRI_hr_ws10 AND minic_ws7_10m AND planc_ws13_hr_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_60cells_hr AND geom_10m_fl1_L17 AND crosc_ws7_10m AND terraintexture_hr_ws49_tp25 AND terraintexture_hr_ws41_tp25 AND maxic_ws13_hr_hr AND planc_DTM_50m_avg_ws11_50m AND vectorruggedness_hr_ws59 AND geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_40cells_hr AND Convexity_10m AND maxic_ws7_10m AND sagaTopographic_Wetness_Index_50m AND geom_hr_L3_fl1_rplipdominance_UE_hr_10cells_hr AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_20cells_hr AND Mass_Balance_Index_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr AND TRI_hr_ws44 AND slope_ws5_hr_hr AND Tangential_Curvature_hr AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_5cells_hr AND Profile_Curvature_50m AND LongitudinalCurvature_10m AND geom_hr_L50m_fl1_rplipmps_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_5cells_hr AND Vertical_Distance_to_Channel_Network_50m"
##
## preds 1 2 3 4 5
## 1 47 9 5 4 0
## 2 0 3 0 0 0
## 3 0 0 10 0 1
## 4 0 2 2 18 2
## 5 0 0 0 0 5
## [1] "Kappa overall = 0.654377880184332"
## [1] "Tau overall = 0.710648148148148"
## [1] "mean quality = 0.552155067155067"
## [1] "The quality is 0.552155067155067"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.686929444529617"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.537037037037037"
##
## preds 1 2 3 4 5
## 1 47 12 12 15 5
## 2 0 0 0 0 0
## 3 0 0 3 0 0
## 4 0 2 2 7 3
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.216277746158224"
## [1] "Tau overall = 0.409722222222222"
## [1] "mean quality = 0.186866683012728"
## [1] "The quality is 0.186866683012728"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.5 - Potential as buffer for acidic contaminants"
## [1] "Prediction error at end is: 0.61038961038961"
## [2] "Prediction error at end is: 0.592207792207792"
## [3] "Prediction error at end is: 0.611688311688312"
## [4] "Prediction error at end is: 0.612554112554113"
## [5] "Prediction error at end is: 0.592640692640693"
## [6] "Prediction error at end is: 0.647619047619048"
## [7] "Prediction error at end is: 0.601298701298701"
## [8] "Prediction error at end is: 0.638961038961039"
## [9] "Prediction error at end is: 0.638528138528139"
## [10] "Prediction error at end is: 0.574891774891775"
## [1] "10fold cv-error: 0.518518518518519"
## [1] "For predictors Standardized_Height_50m AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_10cells_hr"
##
## preds 1 2 3 4 5
## 1 38 12 9 7 3
## 2 0 10 1 2 0
## 3 0 1 10 3 5
## 4 1 1 0 4 1
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.391847227322806"
## [1] "Tau overall = 0.467592592592593"
## [1] "mean quality = 0.293716283044777"
## [1] "The quality is 0.293716283044777"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.537037037037037"
## [1] "For predictors Standardized_Height_50m AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_60cells_hr AND planc_DTM_50m_avg_ws3_50m AND geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr AND Normalized_Height_50m AND Mid_Slope_Positon_50m AND Slope_Height_50m AND LS_Factor_hr AND Profile_Curvature_50m AND ValleyDepth_10m AND Texture_10m AND geom_hr_L3_fl1_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_40cells_hr AND longc_ws11_10m AND Texture_50m AND terraintexture_hr_ws41_tp25 AND fischerk_ws33_hr AND Catchment_slope_50m AND Convexity_50m AND sagaTopographic_Wetness_Index_50m AND geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr AND geom_hr_L3_fl10_rplipshape_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND planc_DTM_50m_avg_ws7_50m AND TRI_hr_ws48 AND dtm_hr_CONVEX_r5 AND General_Curvature_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_5cells_hr AND geom_hr_L3_fl1_rplipshannon_UE_hr_5cells_hr AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_10cells_hr AND longc_DTM_50m_avg_ws7_50m AND geom_hr_L3_fl1_rpliprichness_UE_hr_60cells_hr AND profc_ws23_hr_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_40cells_hr AND planc_ws23_hr_hr AND geom_hr_L50m_fl10_rplipsimpson_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_10cells_hr AND crosc_ws3_hr_hr AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr"
##
## preds 1 2 3 4 5
## 1 38 4 4 0 1
## 2 1 20 2 1 2
## 3 0 0 14 0 3
## 4 0 0 0 15 0
## 5 0 0 0 0 3
## [1] "Kappa overall = 0.773953488372093"
## [1] "Tau overall = 0.791666666666667"
## [1] "mean quality = 0.667572463768116"
## [1] "The quality is 0.667572463768116"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.774060526135238"
## [1] "######################WITH ALL PREDICTORS###################"
## [1] "10fold cv-error: 0.638888888888889"
##
## preds 1 2 3 4 5
## 1 39 21 20 16 9
## 2 0 3 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 0 0
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.0492196878751501"
## [1] "Tau overall = 0.236111111111111"
## [1] "mean quality = 0.0992857142857143"
## [1] "The quality is 0.0992857142857143"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
i=1
for(i in 1:length(dependentlist)){
dep=dependentlist[i]
print(dependentlist_eng[i])
load(paste("../data/modeldata/SVMorigmodeldatawithgeoandgeom_",dep,".RData",sep=""))
preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/SVM_localterrain/SVMwithgeoandgeom_fw_5fold_10p_",dep,"_localterrain",sep=""),kk=1:5,endround = 10,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
predict_radial_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,kappasum = F,tausum = F)
predictors = vector()
for (n in 1:ncol(preds)){
for (n2 in 1:nrow(preds)){
predictors <- c(predictors,as.character(preds[n2,n]))
}
}
print("###################### WITH PREDICTORS from the FW SELECTION ###################")
uniquepredictors <- unique(predictors)
predict_radial_full(modeldata=origmodeldata, dependent=dep, predictors=uniquepredictors)
#print("######################WITH ALL PREDICTORS###################")
#predict_radial_full(modeldata=origmodeldata, dependent=dep, predictors=unlist(paramsets[1]),printpreds = FALSE)
rm(origmodeldata,paramsets,paramsetnames,dependent)
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"
## [1] "Prediction error at end is: 0.657575757575758"
## [2] "Prediction error at end is: 0.686580086580087"
## [3] "Prediction error at end is: 0.658441558441558"
## [4] "Prediction error at end is: 0.648917748917749"
## [5] "Prediction error at end is: 0.667099567099567"
## [6] "Prediction error at end is: 0.694805194805195"
## [7] "Prediction error at end is: 0.694372294372294"
## [8] "Prediction error at end is: 0.703896103896104"
## [9] "Prediction error at end is: 0.666233766233766"
## [10] "Prediction error at end is: 0.676190476190476"
## [1] "10fold cv-error: 0.5"
## [1] "For predictors slope_DTM_50m_avg_ws9_50m AND planc_ws19_hr_hr"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 0 0 0 0 0
## 3 1 0 4 0 1
## 4 3 11 12 28 6
## 5 0 2 5 10 25
## [1] "Kappa overall = 0.304018195602729"
## [1] "Tau overall = 0.409722222222222"
## [1] "mean quality = 0.216823425022183"
## [1] "The quality is 0.216823425022183"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.62962962962963"
## [1] "For predictors slope_DTM_50m_avg_ws9_50m AND planc_ws19_hr_hr AND Longitudinal_Curvature_hr AND profc_DTM_50m_avg_ws11_50m AND maxic_ws15_10m AND GeneralCurvature_10m AND maxic_DTM_50m_avg_ws9_50m AND profc_DTM_50m_avg_ws9_50m AND MaximalCurvature_10m AND slope_ws5_hr_hr AND slope_DTM_50m_avg_ws11_50m AND profc_ws3_10m AND planc_ws11_hr_hr AND crosc_ws7_hr_hr AND crosc_ws3_hr_hr AND slope_DTM_50m_avg_ws3_50m AND ProfileCurvature_10m AND slope_ws5_10m AND longc_ws3_10m AND profc_ws5_hr_hr AND crosc_ws3_10m AND minic_ws19_hr_hr AND minic_ws3_10m AND minic_ws29_hr_hr AND profc_ws7_hr_hr AND minic_ws5_hr_hr AND crosc_ws5_10m AND minic_ws13_hr_hr AND maxic_ws5_hr_hr AND slope_ws3_hr_hr AND TangentialCurvature_10m AND CrossSectionalCurvature_50m AND planc_ws5_10m AND slope_DTM_50m_avg_ws7_50m AND longc_DTM_50m_avg_ws3_50m AND longc_ws7_10m AND Flow_Line_Curvature_hr AND slope_ws19_hr_hr AND profc_ws13_hr_hr"
##
## preds 1 2 3 4 5
## 1 2 0 0 0 0
## 2 0 7 0 0 0
## 3 1 2 9 0 0
## 4 1 4 9 33 5
## 5 0 0 3 5 27
## [1] "Kappa overall = 0.605695509309967"
## [1] "Tau overall = 0.652777777777778"
## [1] "mean quality = 0.533481781376518"
## [1] "The quality is 0.533481781376518"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.662776125620467"
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"
## [1] "Prediction error at end is: 0.667099567099567"
## [2] "Prediction error at end is: 0.648484848484848"
## [3] "Prediction error at end is: 0.638961038961039"
## [4] "Prediction error at end is: 0.61991341991342"
## [5] "Prediction error at end is: 0.619047619047619"
## [6] "Prediction error at end is: 0.573160173160173"
## [7] "Prediction error at end is: 0.583549783549783"
## [8] "Prediction error at end is: 0.545454545454545"
## [9] "Prediction error at end is: 0.535930735930736"
## [10] "Prediction error at end is: 0.564069264069264"
## [1] "10fold cv-error: 0.564814814814815"
## [1] "For predictors crosc_ws5_10m AND CrossSectionalCurvature_hr"
##
## preds 2 3 4 5
## 2 32 20 22 5
## 3 0 2 0 0
## 4 1 1 13 5
## 5 0 0 0 7
## [1] "Kappa overall = 0.2880859375"
## [1] "Tau overall = 0.333333333333333"
## [1] "mean quality = 0.302061259286323"
## [1] "The quality is 0.302061259286323"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.456055648375919"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.574074074074074"
## [1] "For predictors crosc_ws5_10m AND CrossSectionalCurvature_hr AND minic_ws5_hr_hr AND minic_ws29_hr_hr AND longc_DTM_50m_avg_ws11_50m AND planc_DTM_50m_avg_ws5_50m AND slope_ws3_10m AND dtm_hr_CONVEX_r30 AND minic_DTM_50m_avg_ws11_50m AND Longitudinal_Curvature_hr AND slope_DTM_50m_avg_ws7_50m AND maxic_ws11_hr_hr AND maxic_ws13_hr_hr AND maxic_ws9_hr_hr AND maxic_ws5_10m AND maxic_ws23_hr_hr AND maxic_ws3_10m AND Convexity_50m AND planc_ws3_hr_hr AND planc_ws7_10m AND Convergence_Index_hr AND slope_DTM_50m_avg_ws11_50m AND crosc_ws11_10m AND Convergence_Index_50m AND slope_DTM_50m_avg_ws9_50m AND minic_ws7_10m AND planc_ws29_hr_hr AND maxic_ws7_10m AND LongitudinalCurvature_10m AND slope_ws5_hr_hr AND crosc_ws5_hr_hr AND General_Curvature_50m AND planc_ws5_10m AND longc_ws5_hr_hr AND maxic_ws11_10m AND slope_ws11_hr_hr"
##
## preds 2 3 4 5
## 2 30 11 5 2
## 3 0 10 0 0
## 4 3 2 30 4
## 5 0 0 0 11
## [1] "Kappa overall = 0.648590021691974"
## [1] "Tau overall = 0.666666666666667"
## [1] "mean quality = 0.587973727040223"
## [1] "The quality is 0.587973727040223"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.701030983148723"
## [1] "1a.3 - Habitat for soil organisms"
## [1] "Prediction error at end is: 0.638961038961039"
## [2] "Prediction error at end is: 0.602164502164502"
## [3] "Prediction error at end is: 0.566666666666667"
## [4] "Prediction error at end is: 0.576190476190476"
## [5] "Prediction error at end is: 0.520779220779221"
## [6] "Prediction error at end is: 0.52034632034632"
## [7] "Prediction error at end is: 0.539393939393939"
## [8] "Prediction error at end is: 0.511688311688312"
## [9] "Prediction error at end is: 0.52987012987013"
## [10] "Prediction error at end is: 0.52987012987013"
## [1] "10fold cv-error: 0.472222222222222"
## [1] "For predictors slope_DTM_50m_avg_ws9_50m AND crosc_ws3_hr_hr"
##
## preds 1 2 3 4
## 1 3 0 0 0
## 2 21 33 7 10
## 3 0 1 18 4
## 4 0 0 1 10
## [1] "Kappa overall = 0.428571428571428"
## [1] "Tau overall = 0.45679012345679"
## [1] "mean quality = 0.390994623655914"
## [1] "The quality is 0.390994623655914"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.5399090352713"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.5"
## [1] "For predictors slope_DTM_50m_avg_ws9_50m AND crosc_ws3_hr_hr AND minic_ws5_10m AND PlanCurvature_10m AND slope_DTM_50m_avg_ws11_50m AND profc_ws7_10m AND maxic_DTM_50m_avg_ws5_50m AND Total_Curvature_50m AND maxic_ws3_10m AND FlowLineCurvature_10m AND slope_ws5_hr_hr AND longc_ws7_10m AND planc_ws19_hr_hr AND slope_DTM_50m_avg_ws7_50m AND minic_ws5_hr_hr AND Flow_Line_Curvature_hr AND slope_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws7_50m AND minic_DTM_50m_avg_ws11_50m AND crosc_ws11_hr_hr AND planc_ws5_10m AND Longitudinal_Curvature_hr AND Convexity_50m AND minic_DTM_50m_avg_ws7_50m AND longc_ws5_hr_hr AND maxic_DTM_50m_avg_ws3_50m AND longc_DTM_50m_avg_ws5_50m AND planc_ws7_10m AND Slope_50m AND maxic_DTM_50m_avg_ws11_50m AND dtm_hr_CONVEX_r20 AND crosc_ws11_10m AND crosc_DTM_50m_avg_ws11_50m AND minic_ws29_hr_hr AND LongitudinalCurvature_10m AND planc_ws3_hr_hr"
##
## preds 1 2 3 4
## 1 18 7 2 1
## 2 6 25 4 5
## 3 0 2 19 1
## 4 0 0 1 17
## [1] "Kappa overall = 0.636995827538248"
## [1] "Tau overall = 0.641975308641975"
## [1] "mean quality = 0.59369706503291"
## [1] "The quality is 0.59369706503291"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.680333202532746"
## [1] "1a.4 - Habitat for crops"
## [1] "Prediction error at end is: 0.176190476190476"
## [2] "Prediction error at end is: 0.222943722943723"
## [3] "Prediction error at end is: 0.213852813852814"
## [4] "Prediction error at end is: 0.223376623376623"
## [5] "Prediction error at end is: 0.204761904761905"
## [6] "Prediction error at end is: 0.213419913419913"
## [7] "Prediction error at end is: 0.222943722943723"
## [8] "Prediction error at end is: 0.232467532467532"
## [9] "Prediction error at end is: 0.232900432900433"
## [10] "Prediction error at end is: 0.241991341991342"
## [1] "10fold cv-error: 0.175925925925926"
## [1] "For predictors slope_ws3_hr_hr AND planc_ws29_hr_hr"
##
## preds 3 4 5
## 3 2 0 0
## 4 5 61 6
## 5 1 1 32
## [1] "Kappa overall = 0.761710794297352"
## [1] "Tau overall = 0.819444444444444"
## [1] "mean quality = 0.628538812785388"
## [1] "The quality is 0.628538812785388"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.684411620853101"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.194444444444444"
## [1] "For predictors slope_ws3_hr_hr AND planc_ws29_hr_hr AND slope_ws15_10m AND planc_ws23_hr_hr AND Plan_Curvature_hr AND Slope_hr AND planc_ws3_hr_hr AND planc_DTM_50m_avg_ws11_50m AND Convergence_Index_50m AND planc_DTM_50m_avg_ws7_50m AND slope_ws11_10m AND slope_DTM_50m_avg_ws3_50m AND Slope_50m AND Convexity_10m AND longc_DTM_50m_avg_ws11_50m AND crosc_ws5_10m AND crosc_ws23_hr_hr AND Flow_Line_Curvature_50m AND Total_Curvature_hr AND minic_ws5_10m AND slope_ws5_hr_hr AND minic_ws7_10m AND longc_ws7_hr_hr AND minic_ws7_hr_hr AND minic_ws15_10m AND slope_ws3_10m AND slope_ws5_10m AND slope_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws11_50m AND slope_ws11_hr_hr AND dtm_hr_CONVEX_r30 AND slope_ws15_hr_hr AND LongitudinalCurvature_10m AND minic_ws11_10m AND planc_ws9_hr_hr AND planc_ws5_10m AND slope_DTM_50m_avg_ws9_50m AND planc_ws11_10m"
##
## preds 3 4 5
## 3 2 0 0
## 4 6 61 6
## 5 0 1 32
## [1] "Kappa overall = 0.760736196319018"
## [1] "Tau overall = 0.819444444444444"
## [1] "mean quality = 0.631612381612382"
## [1] "The quality is 0.631612381612382"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.694203111304538"
## [1] "1c.1 - Average precipitation retention capacity"
## [1] "Prediction error at end is: 0.611255411255411"
## [2] "Prediction error at end is: 0.583982683982684"
## [3] "Prediction error at end is: 0.583982683982684"
## [4] "Prediction error at end is: 0.556709956709957"
## [5] "Prediction error at end is: 0.592640692640693"
## [6] "Prediction error at end is: 0.62034632034632"
## [7] "Prediction error at end is: 0.611255411255411"
## [8] "Prediction error at end is: 0.611255411255411"
## [9] "Prediction error at end is: 0.611255411255411"
## [10] "Prediction error at end is: 0.611255411255411"
## [1] "10fold cv-error: 0.518518518518519"
## [1] "For predictors slope_DTM_50m_avg_ws5_50m AND Flow_Line_Curvature_50m"
##
## preds 1 2 3 4 5
## 1 11 1 3 3 5
## 2 15 44 9 14 2
## 3 0 0 0 0 0
## 4 0 0 0 1 0
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.227297743533297"
## [1] "Tau overall = 0.398148148148148"
## [1] "mean quality = 0.172535259717922"
## [1] "The quality is 0.172535259717922"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.537037037037037"
## [1] "For predictors slope_DTM_50m_avg_ws5_50m AND Flow_Line_Curvature_50m AND maxic_DTM_50m_avg_ws11_50m AND dtm_hr_CONVEX_r10 AND crosc_ws3_10m AND slope_ws11_10m AND slope_DTM_50m_avg_ws7_50m AND Plan_Curvature_hr AND planc_ws13_hr_hr AND slope_ws15_10m AND ProfileCurvature_10m AND slope_DTM_50m_avg_ws11_50m AND crosc_DTM_50m_avg_ws11_50m AND crosc_DTM_50m_avg_ws9_50m AND planc_ws5_hr_hr AND dtm_hr_CONVEX_r30 AND planc_ws9_hr_hr AND minic_ws11_hr_hr AND profc_ws15_10m AND minic_DTM_50m_avg_ws7_50m AND minic_ws19_hr_hr AND minic_DTM_50m_avg_ws9_50m AND minic_ws3_hr_hr AND minic_DTM_50m_avg_ws11_50m AND profc_DTM_50m_avg_ws9_50m AND profc_ws5_hr_hr AND CrossSectionalCurvature_50m AND crosc_ws19_hr_hr AND GeneralCurvature_10m AND longc_ws11_10m AND profc_ws3_10m AND minic_ws15_10m AND slope_ws7_10m AND minic_DTM_50m_avg_ws5_50m AND slope_ws29_hr_hr AND Total_Curvature_50m AND planc_ws15_hr_hr AND Total_Curvature_hr AND planc_DTM_50m_avg_ws5_50m AND planc_ws11_hr_hr AND maxic_DTM_50m_avg_ws3_50m AND slope_DTM_50m_avg_ws9_50m"
##
## preds 1 2 3 4 5
## 1 16 0 3 1 3
## 2 10 45 8 12 2
## 3 0 0 1 0 0
## 4 0 0 0 5 0
## 5 0 0 0 0 2
## [1] "Kappa overall = 0.437274549098196"
## [1] "Tau overall = 0.548611111111111"
## [1] "mean quality = 0.343217893217893"
## [1] "The quality is 0.343217893217893"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.495174407613578"
## [1] "1c.1 - Minimum precipitation retention capacity"
## [1] "Prediction error at end is: 0.75974025974026"
## [2] "Prediction error at end is: 0.767532467532468"
## [3] "Prediction error at end is: 0.767099567099567"
## [4] "Prediction error at end is: 0.738961038961039"
## [5] "Prediction error at end is: 0.73030303030303"
## [6] "Prediction error at end is: 0.738961038961039"
## [7] "Prediction error at end is: 0.738528138528138"
## [8] "Prediction error at end is: 0.766666666666667"
## [9] "Prediction error at end is: 0.767099567099567"
## [10] "Prediction error at end is: 0.767099567099567"
## [1] "10fold cv-error: 0.703703703703704"
## [1] "For predictors crosc_ws11_10m AND minic_ws15_hr_hr"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 6 21 16 15 6
## 3 0 0 1 0 0
## 4 4 4 4 10 2
## 5 2 1 1 1 14
## [1] "Kappa overall = 0.250671441360788"
## [1] "Tau overall = 0.282407407407407"
## [1] "mean quality = 0.223664178012004"
## [1] "The quality is 0.223664178012004"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.731481481481482"
## [1] "For predictors crosc_ws11_10m AND minic_ws15_hr_hr AND planc_ws5_10m AND minic_DTM_50m_avg_ws11_50m AND Total_Curvature_50m AND Plan_Curvature_50m AND crosc_DTM_50m_avg_ws3_50m AND planc_ws3_hr_hr AND planc_ws11_10m AND Maximal_Curvature_50m AND minic_ws11_10m AND profc_ws15_10m AND Longitudinal_Curvature_hr AND PlanCurvature_10m AND Profile_Curvature_50m AND Convergence_Index_hr AND minic_DTM_50m_avg_ws5_50m AND profc_DTM_50m_avg_ws11_50m AND planc_ws7_10m AND CrossSectionalCurvature_50m AND minic_DTM_50m_avg_ws3_50m AND crosc_ws15_hr_hr AND dtm_hr_CONVEX_r30 AND slope_DTM_50m_avg_ws9_50m AND planc_ws3_10m AND maxic_ws29_hr_hr AND maxic_ws5_10m AND maxic_ws7_hr_hr AND maxic_ws3_10m AND longc_ws29_hr_hr AND Flow_Line_Curvature_50m AND minic_ws13_hr_hr AND slope_ws11_10m AND minic_DTM_50m_avg_ws7_50m AND LongitudinalCurvature_10m AND MinimalCurvature_10m AND longc_DTM_50m_avg_ws3_50m AND Profile_Curvature_hr AND Convexity_10m AND profc_ws11_10m AND minic_ws5_hr_hr AND planc_ws29_hr_hr"
##
## preds 1 2 3 4 5
## 1 3 0 0 0 0
## 2 4 25 8 7 3
## 3 1 1 12 2 1
## 4 1 0 1 16 3
## 5 3 0 1 1 15
## [1] "Kappa overall = 0.558257793499889"
## [1] "Tau overall = 0.571759259259259"
## [1] "mean quality = 0.45739247311828"
## [1] "The quality is 0.45739247311828"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.581372327644911"
## [1] "1c.2 - Retention capacity for heavy precipitation events"
## [1] "Prediction error at end is: 0.333333333333333"
## [2] "Prediction error at end is: 0.296536796536797"
## [3] "Prediction error at end is: 0.287445887445887"
## [4] "Prediction error at end is: 0.287445887445887"
## [5] "Prediction error at end is: 0.287445887445887"
## [6] "Prediction error at end is: 0.278354978354978"
## [7] "Prediction error at end is: 0.278354978354978"
## [8] "Prediction error at end is: 0.287445887445887"
## [9] "Prediction error at end is: 0.287445887445887"
## [10] "Prediction error at end is: 0.278354978354978"
## [1] "10fold cv-error: 0.268518518518518"
## [1] "For predictors crosc_ws7_10m AND maxic_ws11_10m"
##
## preds 1 2 3 4 5
## 1 75 10 6 3 8
## 2 0 1 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 5 0
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.264193792581378"
## [1] "Tau overall = 0.6875"
## [1] "mean quality = 0.29024064171123"
## [1] "The quality is 0.29024064171123"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.287037037037037"
## [1] "For predictors crosc_ws7_10m AND maxic_ws11_10m AND Total_Curvature_50m AND slope_DTM_50m_avg_ws3_50m AND crosc_DTM_50m_avg_ws5_50m AND crosc_ws3_hr_hr AND slope_ws11_10m AND Profile_Curvature_50m AND longc_DTM_50m_avg_ws3_50m AND profc_ws19_hr_hr AND longc_ws15_hr_hr AND Minimal_Curvature_50m AND Minimal_Curvature_hr AND minic_ws15_hr_hr AND profc_ws15_hr_hr AND longc_ws13_hr_hr AND longc_ws11_hr_hr AND minic_ws19_hr_hr AND planc_ws11_10m AND crosc_ws11_10m AND Maximal_Curvature_50m AND Convexity_50m AND Flow_Line_Curvature_50m AND Total_Curvature_hr AND Plan_Curvature_hr AND slope_DTM_50m_avg_ws5_50m AND profc_ws11_10m AND longc_ws5_10m AND minic_ws3_hr_hr AND maxic_ws7_10m AND maxic_ws15_10m AND maxic_DTM_50m_avg_ws5_50m AND Longitudinal_Curvature_hr AND minic_ws29_hr_hr AND minic_ws15_10m AND minic_ws5_10m AND longc_ws15_10m AND maxic_ws23_hr_hr AND minic_ws11_10m AND General_Curvature_hr AND maxic_ws3_hr_hr AND General_Curvature_50m"
##
## preds 1 2 3 4 5
## 1 75 6 6 0 6
## 2 0 5 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 8 0
## 5 0 0 0 0 2
## [1] "Kappa overall = 0.573122529644269"
## [1] "Tau overall = 0.791666666666667"
## [1] "mean quality = 0.502199413489736"
## [1] "The quality is 0.502199413489736"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1c.3 - groundwater reformation rate"
## [1] "Prediction error at end is: 0.658441558441558"
## [2] "Prediction error at end is: 0.649350649350649"
## [3] "Prediction error at end is: 0.658874458874459"
## [4] "Prediction error at end is: 0.630735930735931"
## [5] "Prediction error at end is: 0.639393939393939"
## [6] "Prediction error at end is: 0.658441558441558"
## [7] "Prediction error at end is: 0.658441558441558"
## [8] "Prediction error at end is: 0.658441558441558"
## [9] "Prediction error at end is: 0.658441558441558"
## [10] "Prediction error at end is: 0.649350649350649"
## [1] "10fold cv-error: 0.537037037037037"
## [1] "For predictors minic_ws15_10m AND maxic_ws7_hr_hr"
##
## preds 1 2 3 4 5
## 1 1 0 0 0 0
## 2 13 37 10 14 9
## 3 0 0 1 0 0
## 4 0 0 1 4 0
## 5 0 2 2 1 13
## [1] "Kappa overall = 0.289833080424886"
## [1] "Tau overall = 0.398148148148148"
## [1] "mean quality = 0.251926548397137"
## [1] "The quality is 0.251926548397137"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.38726358510718"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.62962962962963"
## [1] "For predictors minic_ws15_10m AND maxic_ws7_hr_hr AND minic_ws29_hr_hr AND minic_ws3_hr_hr AND profc_ws3_hr_hr AND minic_ws7_10m AND LongitudinalCurvature_10m AND planc_DTM_50m_avg_ws11_50m AND crosc_ws5_10m AND maxic_ws3_10m AND slope_DTM_50m_avg_ws7_50m AND slope_ws3_10m AND Convexity_50m AND FlowLineCurvature_10m AND slope_DTM_50m_avg_ws5_50m AND crosc_ws7_10m AND dtm_hr_CONVEX_r10 AND General_Curvature_hr AND slope_ws7_10m AND crosc_ws23_hr_hr AND CrossSectionalCurvature_hr AND crosc_ws19_hr_hr AND planc_ws29_hr_hr AND planc_DTM_50m_avg_ws3_50m AND Flow_Line_Curvature_hr AND maxic_ws3_hr_hr AND crosc_ws3_hr_hr AND slope_DTM_50m_avg_ws9_50m AND Tangential_Curvature_hr AND slope_ws15_10m AND Minimal_Curvature_hr AND minic_ws11_10m AND minic_DTM_50m_avg_ws3_50m AND profc_ws7_hr_hr AND longc_DTM_50m_avg_ws9_50m AND Profile_Curvature_hr AND minic_ws13_hr_hr"
##
## preds 1 2 3 4 5
## 1 3 0 0 0 0
## 2 11 38 10 5 9
## 3 0 0 2 0 0
## 4 0 0 1 12 0
## 5 0 1 1 2 13
## [1] "Kappa overall = 0.468373123307901"
## [1] "Tau overall = 0.537037037037037"
## [1] "mean quality = 0.394131274131274"
## [1] "The quality is 0.394131274131274"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.549432082090476"
## [1] "1c.4 - Potential for providing nutrients for plants"
## [1] "Prediction error at end is: 0.352380952380952"
## [2] "Prediction error at end is: 0.388311688311688"
## [3] "Prediction error at end is: 0.36969696969697"
## [4] "Prediction error at end is: 0.341991341991342"
## [5] "Prediction error at end is: 0.351948051948052"
## [6] "Prediction error at end is: 0.370995670995671"
## [7] "Prediction error at end is: 0.352380952380952"
## [8] "Prediction error at end is: 0.343290043290043"
## [9] "Prediction error at end is: 0.343290043290043"
## [10] "Prediction error at end is: 0.343290043290043"
## [1] "10fold cv-error: 0.342592592592593"
## [1] "For predictors minic_ws15_hr_hr AND crosc_ws3_hr_hr"
##
## preds 1 3 5
## 1 64 9 15
## 3 0 0 0
## 5 1 2 17
## [1] "Kappa overall = 0.450226244343891"
## [1] "Tau overall = 0.625"
## [1] "mean quality = 0.401605136436597"
## [1] "The quality is 0.401605136436597"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.287037037037037"
## [1] "For predictors minic_ws15_hr_hr AND crosc_ws3_hr_hr AND minic_ws5_hr_hr AND maxic_ws3_10m AND planc_DTM_50m_avg_ws7_50m AND planc_DTM_50m_avg_ws11_50m AND Slope_50m AND minic_ws7_hr_hr AND Convergence_Index_50m AND Convergence_Index_hr AND minic_ws9_hr_hr AND Profile_Curvature_50m AND profc_ws23_hr_hr AND ProfileCurvature_10m AND Convexity_10m AND Total_Curvature_hr AND profc_ws15_hr_hr AND Tangential_Curvature_50m AND dtm_hr_CONVEX_r30 AND slope_ws3_hr_hr AND slope_ws15_hr_hr AND crosc_ws3_10m AND Plan_Curvature_50m AND crosc_DTM_50m_avg_ws5_50m AND General_Curvature_50m AND longc_ws15_10m AND minic_DTM_50m_avg_ws11_50m AND planc_ws11_10m AND Total_Curvature_50m AND Minimal_Curvature_50m AND LongitudinalCurvature_10m AND minic_ws7_10m AND General_Curvature_hr AND planc_ws3_hr_hr AND slope_ws5_hr_hr AND profc_ws29_hr_hr AND minic_ws15_10m AND Flow_Line_Curvature_hr"
##
## preds 1 3 5
## 1 64 7 6
## 3 0 2 0
## 5 1 2 26
## [1] "Kappa overall = 0.697320021019443"
## [1] "Tau overall = 0.777777777777778"
## [1] "mean quality = 0.581729381729382"
## [1] "The quality is 0.581729381729382"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.638035811797989"
## [1] "1c.5 - Potential as a CO2 sink"
## [1] "Prediction error at end is: 0.502164502164502"
## [2] "Prediction error at end is: 0.446320346320346"
## [3] "Prediction error at end is: 0.38008658008658"
## [4] "Prediction error at end is: 0.351948051948052"
## [5] "Prediction error at end is: 0.37012987012987"
## [6] "Prediction error at end is: 0.37965367965368"
## [7] "Prediction error at end is: 0.370995670995671"
## [8] "Prediction error at end is: 0.370995670995671"
## [9] "Prediction error at end is: 0.361904761904762"
## [10] "Prediction error at end is: 0.370995670995671"
## [1] "10fold cv-error: 0.333333333333333"
## [1] "For predictors minic_ws11_hr_hr AND maxic_ws3_10m"
##
## preds 1 2 3 4 5
## 1 44 0 3 2 1
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 7 1 15 31 4
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.505"
## [1] "Tau overall = 0.618055555555555"
## [1] "mean quality = 0.257719298245614"
## [1] "The quality is 0.257719298245614"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.314814814814815"
## [1] "For predictors minic_ws11_hr_hr AND maxic_ws3_10m AND MaximalCurvature_10m AND maxic_ws3_hr_hr AND maxic_ws19_hr_hr AND Convexity_10m AND dtm_hr_CONVEX_r30 AND slope_ws3_10m AND Minimal_Curvature_50m AND planc_DTM_50m_avg_ws7_50m AND crosc_ws13_hr_hr AND minic_ws13_hr_hr AND minic_ws5_hr_hr AND ProfileCurvature_10m AND crosc_DTM_50m_avg_ws5_50m AND planc_ws3_hr_hr AND crosc_ws5_hr_hr AND maxic_ws5_10m AND minic_DTM_50m_avg_ws5_50m AND profc_ws3_10m AND longc_ws3_hr_hr AND minic_ws15_10m AND minic_ws9_hr_hr AND longc_ws15_10m AND Maximal_Curvature_hr AND Profile_Curvature_50m AND slope_ws15_hr_hr AND planc_ws5_hr_hr AND slope_DTM_50m_avg_ws7_50m AND crosc_ws15_hr_hr AND profc_ws5_hr_hr AND planc_DTM_50m_avg_ws11_50m AND longc_ws3_10m AND Slope_hr AND TangentialCurvature_10m AND minic_ws3_hr_hr AND minic_ws7_hr_hr AND slope_DTM_50m_avg_ws11_50m AND planc_DTM_50m_avg_ws9_50m AND longc_DTM_50m_avg_ws9_50m"
##
## preds 1 2 3 4 5
## 1 46 0 1 1 1
## 2 0 0 0 0 0
## 3 0 0 6 0 0
## 4 5 1 11 32 4
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.645320197044335"
## [1] "Tau overall = 0.722222222222222"
## [1] "mean quality = 0.355555555555556"
## [1] "The quality is 0.355555555555556"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.1 - Potential for retention of heavy metals"
## [1] "Prediction error at end is: 0.425974025974026"
## [2] "Prediction error at end is: 0.425541125541126"
## [3] "Prediction error at end is: 0.388744588744589"
## [4] "Prediction error at end is: 0.369264069264069"
## [5] "Prediction error at end is: 0.379220779220779"
## [6] "Prediction error at end is: 0.37012987012987"
## [7] "Prediction error at end is: 0.351515151515152"
## [8] "Prediction error at end is: 0.37012987012987"
## [9] "Prediction error at end is: 0.37965367965368"
## [10] "Prediction error at end is: 0.37012987012987"
## [1] "10fold cv-error: 0.37962962962963"
## [1] "For predictors crosc_ws13_hr_hr AND slope_ws19_hr_hr"
##
## preds 1 2 3 4 5
## 1 38 5 8 3 11
## 2 0 0 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 0 0
## 5 4 2 2 4 31
## [1] "Kappa overall = 0.409090909090909"
## [1] "Tau overall = 0.548611111111111"
## [1] "mean quality = 0.224959742351047"
## [1] "The quality is 0.224959742351047"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.342592592592593"
## [1] "For predictors crosc_ws13_hr_hr AND slope_ws19_hr_hr AND minic_ws11_10m AND maxic_ws3_10m AND planc_ws3_hr_hr AND minic_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND slope_ws7_10m AND slope_DTM_50m_avg_ws3_50m AND maxic_ws11_hr_hr AND slope_ws11_10m AND dtm_hr_CONVEX_r30 AND longc_DTM_50m_avg_ws9_50m AND TotalCurvature_10m AND Total_Curvature_50m AND minic_ws23_hr_hr AND longc_DTM_50m_avg_ws5_50m AND maxic_ws5_hr_hr AND minic_ws9_hr_hr AND minic_ws5_hr_hr AND crosc_DTM_50m_avg_ws5_50m AND crosc_ws3_hr_hr AND maxic_ws3_hr_hr AND slope_DTM_50m_avg_ws5_50m AND maxic_ws15_hr_hr AND slope_ws5_10m AND maxic_ws23_hr_hr AND minic_ws29_hr_hr AND longc_ws9_hr_hr AND Convexity_10m AND General_Curvature_50m AND planc_DTM_50m_avg_ws11_50m AND Slope_10m AND longc_ws19_hr_hr AND Minimal_Curvature_50m AND profc_DTM_50m_avg_ws9_50m AND maxic_ws7_10m AND profc_DTM_50m_avg_ws5_50m AND slope_DTM_50m_avg_ws11_50m"
##
## preds 1 2 3 4 5
## 1 41 4 7 3 5
## 2 0 1 0 0 0
## 3 0 0 2 0 0
## 4 0 0 0 2 0
## 5 1 2 1 2 37
## [1] "Kappa overall = 0.629984925311772"
## [1] "Tau overall = 0.710648148148148"
## [1] "mean quality = 0.414307181889149"
## [1] "The quality is 0.414307181889149"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.55548001516698"
## [1] "1d.2 - Potential for transforming organic contaminants"
## [1] "Prediction error at end is: 0.683549783549784"
## [2] "Prediction error at end is: 0.674891774891775"
## [3] "Prediction error at end is: 0.627272727272727"
## [4] "Prediction error at end is: 0.645887445887446"
## [5] "Prediction error at end is: 0.646320346320346"
## [6] "Prediction error at end is: 0.638528138528139"
## [7] "Prediction error at end is: 0.638961038961039"
## [8] "Prediction error at end is: 0.61991341991342"
## [9] "Prediction error at end is: 0.628138528138528"
## [10] "Prediction error at end is: 0.656277056277056"
## [1] "10fold cv-error: 0.555555555555556"
## [1] "For predictors slope_ws7_10m AND crosc_ws15_10m"
##
## preds 1 2 3 4 5
## 1 0 0 0 0 0
## 2 0 0 0 0 0
## 3 10 18 22 2 9
## 4 0 0 0 3 1
## 5 0 4 6 4 29
## [1] "Kappa overall = 0.292490598083222"
## [1] "Tau overall = 0.375"
## [1] "mean quality = 0.235105604055196"
## [1] "The quality is 0.235105604055196"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.583333333333333"
## [1] "For predictors slope_ws7_10m AND crosc_ws15_10m AND planc_ws5_hr_hr AND Total_Curvature_50m AND minic_ws11_10m AND longc_ws3_hr_hr AND TangentialCurvature_10m AND Slope_hr AND Plan_Curvature_hr AND Convexity_10m AND MaximalCurvature_10m AND profc_DTM_50m_avg_ws7_50m AND planc_DTM_50m_avg_ws9_50m AND profc_DTM_50m_avg_ws5_50m AND maxic_DTM_50m_avg_ws7_50m AND longc_DTM_50m_avg_ws7_50m AND maxic_ws13_hr_hr AND Minimal_Curvature_50m AND minic_ws11_hr_hr AND profc_ws11_hr_hr AND longc_DTM_50m_avg_ws9_50m AND profc_DTM_50m_avg_ws9_50m AND dtm_hr_CONVEX_r30 AND slope_ws29_hr_hr AND dtm_hr_CONVEX_r20 AND maxic_DTM_50m_avg_ws5_50m AND Tangential_Curvature_50m AND minic_ws5_hr_hr AND slope_DTM_50m_avg_ws5_50m AND crosc_ws7_10m AND Slope_10m AND minic_ws19_hr_hr AND slope_DTM_50m_avg_ws7_50m AND Convexity_50m AND dtm_hr_CONVEX_r10 AND minic_DTM_50m_avg_ws3_50m AND crosc_ws11_10m AND slope_DTM_50m_avg_ws9_50m AND minic_ws3_hr_hr AND longc_ws15_hr_hr AND profc_ws3_10m AND Minimal_Curvature_hr AND CrossSectionalCurvature_hr AND planc_ws15_hr_hr"
##
## preds 1 2 3 4 5
## 1 1 0 0 0 0
## 2 0 6 0 0 0
## 3 9 14 27 2 5
## 4 0 0 0 4 0
## 5 0 2 1 3 34
## [1] "Kappa overall = 0.533253301320528"
## [1] "Tau overall = 0.583333333333333"
## [1] "mean quality = 0.407648902821317"
## [1] "The quality is 0.407648902821317"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.583706053095214"
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"
## [1] "Prediction error at end is: 0.0372294372294372"
## [2] "Prediction error at end is: 0.0372294372294372"
## [3] "Prediction error at end is: 0.0372294372294372"
## [4] "Prediction error at end is: 0.0372294372294372"
## [5] "Prediction error at end is: 0.0372294372294372"
## [6] "Prediction error at end is: 0.0372294372294372"
## [7] "Prediction error at end is: 0.0372294372294372"
## [8] "Prediction error at end is: 0.0372294372294372"
## [9] "Prediction error at end is: 0.0372294372294372"
## [10] "Prediction error at end is: 0.0372294372294372"
## [1] "10fold cv-error: 0.0185185185185185"
## [1] "For predictors PlanCurvature_10m AND planc_ws11_hr_hr"
##
## preds 4 5
## 4 2 0
## 5 2 104
## [1] "Kappa overall = 0.658227848101266"
## [1] "Tau overall = 0.962962962962963"
## [1] "mean quality = 0.740566037735849"
## [1] "The quality is 0.740566037735849"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.518568491248852"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.0370370370370371"
## [1] "For predictors PlanCurvature_10m AND planc_ws11_hr_hr AND planc_ws13_hr_hr AND planc_ws19_hr_hr AND Convexity_10m AND planc_ws9_hr_hr AND crosc_DTM_50m_avg_ws11_50m AND maxic_DTM_50m_avg_ws11_50m AND minic_DTM_50m_avg_ws11_50m AND profc_DTM_50m_avg_ws9_50m AND longc_ws7_10m AND profc_ws5_10m AND profc_ws7_10m AND LongitudinalCurvature_10m AND maxic_DTM_50m_avg_ws3_50m AND crosc_ws11_10m AND crosc_ws15_10m AND crosc_ws3_10m AND crosc_ws5_10m AND crosc_ws7_10m AND longc_ws11_10m AND longc_ws15_10m AND longc_ws3_10m AND longc_ws5_10m AND minic_ws15_10m AND minic_ws11_10m AND minic_ws7_10m AND slope_ws15_10m AND Tangential_Curvature_50m"
##
## preds 4 5
## 4 0 0
## 5 4 104
## [1] "Kappa overall = 0"
## [1] "Tau overall = 0.925925925925926"
## [1] "mean quality = 0.481481481481482"
## [1] "The quality is 0.481481481481482"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "1d.4 - Potential for retention of water-soluble contaminants"
## [1] "Prediction error at end is: 0.593073593073593"
## [2] "Prediction error at end is: 0.556709956709957"
## [3] "Prediction error at end is: 0.538095238095238"
## [4] "Prediction error at end is: 0.547186147186147"
## [5] "Prediction error at end is: 0.538095238095238"
## [6] "Prediction error at end is: 0.51991341991342"
## [7] "Prediction error at end is: 0.510822510822511"
## [8] "Prediction error at end is: 0.51991341991342"
## [9] "Prediction error at end is: 0.501298701298701"
## [10] "Prediction error at end is: 0.510822510822511"
## [1] "10fold cv-error: 0.472222222222222"
## [1] "For predictors minic_ws11_10m AND slope_ws5_hr_hr"
##
## preds 1 2 3 4 5
## 1 43 8 7 9 3
## 2 0 0 0 0 0
## 3 0 0 4 0 0
## 4 4 6 6 13 2
## 5 0 0 0 0 3
## [1] "Kappa overall = 0.360526315789474"
## [1] "Tau overall = 0.479166666666667"
## [1] "mean quality = 0.303275039745628"
## [1] "The quality is 0.303275039745628"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.453703703703704"
## [1] "For predictors minic_ws11_10m AND slope_ws5_hr_hr AND Flow_Line_Curvature_50m AND Slope_10m AND Maximal_Curvature_hr AND LongitudinalCurvature_10m AND minic_DTM_50m_avg_ws3_50m AND General_Curvature_hr AND Total_Curvature_50m AND slope_ws5_10m AND maxic_ws29_hr_hr AND slope_ws3_hr_hr AND crosc_ws5_10m AND slope_DTM_50m_avg_ws5_50m AND Flow_Line_Curvature_hr AND minic_ws7_10m AND MaximalCurvature_10m AND General_Curvature_50m AND crosc_ws5_hr_hr AND crosc_ws15_hr_hr AND slope_ws15_10m AND profc_ws5_10m AND planc_ws5_hr_hr AND minic_DTM_50m_avg_ws5_50m AND Total_Curvature_hr AND maxic_ws11_10m AND maxic_DTM_50m_avg_ws11_50m AND slope_ws7_10m AND crosc_ws19_hr_hr AND minic_ws5_hr_hr AND slope_DTM_50m_avg_ws3_50m AND planc_DTM_50m_avg_ws11_50m AND planc_ws9_hr_hr AND longc_ws15_10m AND Slope_hr AND ProfileCurvature_10m AND slope_ws11_hr_hr AND maxic_ws5_10m AND longc_ws3_hr_hr"
##
## preds 1 2 3 4 5
## 1 46 9 6 7 1
## 2 0 2 0 0 0
## 3 0 0 8 0 0
## 4 1 3 3 15 2
## 5 0 0 0 0 5
## [1] "Kappa overall = 0.550526726492392"
## [1] "Tau overall = 0.62962962962963"
## [1] "mean quality = 0.475891840607211"
## [1] "The quality is 0.475891840607211"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.619729938181675"
## [1] "1d.5 - Potential as buffer for acidic contaminants"
## [1] "Prediction error at end is: 0.63030303030303"
## [2] "Prediction error at end is: 0.665800865800866"
## [3] "Prediction error at end is: 0.657142857142857"
## [4] "Prediction error at end is: 0.703463203463203"
## [5] "Prediction error at end is: 0.694372294372294"
## [6] "Prediction error at end is: 0.685714285714286"
## [7] "Prediction error at end is: 0.666233766233766"
## [8] "Prediction error at end is: 0.666666666666667"
## [9] "Prediction error at end is: 0.675757575757576"
## [10] "Prediction error at end is: 0.674891774891775"
## [1] "10fold cv-error: 0.592592592592593"
## [1] "For predictors planc_ws11_10m AND slope_DTM_50m_avg_ws11_50m"
##
## preds 1 2 3 4 5
## 1 39 11 16 12 7
## 2 0 13 2 1 0
## 3 0 0 2 0 1
## 4 0 0 0 3 1
## 5 0 0 0 0 0
## [1] "Kappa overall = 0.297538579262849"
## [1] "Tau overall = 0.409722222222222"
## [1] "mean quality = 0.242402738873327"
## [1] "The quality is 0.242402738873327"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "10fold cv-error: 0.537037037037037"
## [1] "For predictors planc_ws11_10m AND slope_DTM_50m_avg_ws11_50m AND minic_DTM_50m_avg_ws5_50m AND maxic_DTM_50m_avg_ws11_50m AND slope_ws7_10m AND planc_ws7_10m AND profc_ws19_hr_hr AND minic_ws7_hr_hr AND Convexity_10m AND crosc_DTM_50m_avg_ws5_50m AND Profile_Curvature_50m AND planc_ws29_hr_hr AND minic_ws23_hr_hr AND Convexity_50m AND maxic_DTM_50m_avg_ws7_50m AND dtm_hr_CONVEX_r10 AND planc_ws23_hr_hr AND minic_ws5_10m AND planc_ws3_hr_hr AND crosc_DTM_50m_avg_ws9_50m AND planc_ws5_10m AND slope_ws29_hr_hr AND profc_DTM_50m_avg_ws11_50m AND crosc_ws11_10m AND minic_ws3_hr_hr AND slope_DTM_50m_avg_ws3_50m AND planc_ws3_10m AND General_Curvature_50m AND Slope_hr AND profc_DTM_50m_avg_ws3_50m AND Minimal_Curvature_hr AND profc_ws29_hr_hr AND minic_ws15_hr_hr AND profc_ws5_10m AND maxic_DTM_50m_avg_ws5_50m AND crosc_DTM_50m_avg_ws7_50m AND longc_ws23_hr_hr AND minic_DTM_50m_avg_ws3_50m AND longc_ws15_10m AND Plan_Curvature_hr AND Convergence_Index_50m"
##
## preds 1 2 3 4 5
## 1 38 7 11 3 3
## 2 1 17 0 1 0
## 3 0 0 9 0 2
## 4 0 0 0 12 0
## 5 0 0 0 0 4
## [1] "Kappa overall = 0.637497003116759"
## [1] "Tau overall = 0.675925925925926"
## [1] "mean quality = 0.572111222111222"
## [1] "The quality is 0.572111222111222"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.706105756645892"